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ABSTRACT 

The origin and evolution of the X-ray emission in very young stellar objects (YSOs) are 
not yet well understood because it is very hard to observe YSOs in the protostellar phase. We 
study the X-ray properties of Class O-I objects in the Orion Nebula Cluster (ONC) and compare 
them with those of the more evolved Class II and III members. Using Chandra Orion Ultradeep 
Project (COUP) data, we study the X-ray properties of stars in different evolutionary classes: 
luminosities, hydrogen column densities Nh, average plasma temperatures and time variability 
are compared in order to understand if the interaction between the circumstellar material and 
the central object can influence the X-ray emission. We have assembled the deepest and most 
complete photometric catalog of objects in the ONC region from the UV to 8 /im using data 
from the HST Treasury Program, deep and almost simultaneous UBVI and JHK images taken, 
respectively, with WFI@2.2m ESO and ISPI@4m CTIO telescopes, and Spitzer IRAC imaging. 
We select high probability candidate Class O-I protostars, distinguishing between those having 
a spectral energy distribution which rises from K up to 8 /xm (Class 0-Ia) from those where the 
SED rises from K up to 4.5 /im and decreasing afterwards (Class 0-Ib). In addition, we select a 
sample of "bona fide" Class II stars and a set of Class III stars with IR emission consistent with 
normal photospheres. 

Our principal result is that Class 0-Ia objects are significantly less luminous in X-rays, both 
in the total and hard bands, than the more evolved Class II stars with mass larger than 0.5 M©; 
these latter show X-ray luminosities similar to those of Class 0- lb stars. This result supports the 
hypothesis that the onset of X-ray e mission occurs at a ve ry early stage of star formation and is 
in agreement with the result found in Giardino et al. ( 2007t ). Spectral properties of Class O-I stars 
are similar to those of the more evolved Class II and III objects, except for a larger absorption 
likely due to gas in the envelope or disk of the protostellar objects. Our data suggest that the 
three different classes have similar X-ray temporal variability. 



Subject headings: open clusters: general — open clusters: individual (Orion Nebula Cluster) 
stars — stars: pre-main sequence — stars: activity 

1. INTRODUCTION 
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Young stellar objects (YSOs) exhibit an X-ray 
luminosity significantly higher than similar stars 
in the main sequence phase. However, it is still 
unclear at what stage of their formation process 
YSOs begin to emit X-ra ys and how this X-ra 
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Feigelson et al.ll200lD . 

Evidence has accumulated over the past few 
years that among the YSOs, the X-ray lumi- 
nosity function (XLF) of Classical T Tauri stars 
(CTTS) is different from that of Weak T Tauri 
stars (WTTS). Such a difference has been inves- 
tigated by some because the accretion and out- 
flow processes, characteristic of the CTTS phase, 
involve the interaction of ionized material with 
the star and disk magnetic field, and these me- 
dia play a fundamental role in the coronal X-ray 
emission of young stars. Using the Chandra Orion 



Ultrad e ep Project (COUP) d ata (jGetman et al 



2005bh . IPreibisch etaP (|2005l ) find that the ac 
creting stars, as discriminated using the Call 
8542A IR triplet line, are less X-ray luminous than 
non-accretors, while no differences are found for 
stars with {K — L) IR excesses, indicative of the 
presence of a (possibly passive) disk, with respect 
to stars without IR excesses. These results are 
in agreement with th ose found in the previo us 
work in the ONC by iFlaccomio et al.l (|2003bl H) 
who found a strong difference in the X-ray lumi- 
nosity of accreting and non-accreting stars. 

In order to understand the evolution of the X- 
ray emission during the sta r formation process, 
Preibisch fc FeigelsonI (|2005r ) looked for a possi- 
ble relation between X-ray emission and stellar 
age in the ONC. They found, in the age range 
~ 0.1 — lOMyr, a slow decay in the X-ray lumi- 
nosity, Lx. Because, as the star contracts, the 
bolometric luminosity, Lboh decreases faster than 
Lx, on average Lx/Lboi and the X-ray surface flux. 
Fx, increase during this age range. 

In this work, we want to investigate the on- 
set and the evolution of X-ray emission during 
the early phase of the star formation by statisti- 
cally comparing the X-ray properties of the stars 
in different evolutionary phases rather than con- 
sidering individual stellar ages, usually derived 
from theoretical isochrones. We hope in this 
way to avoid some of the problems introduced 
when using inferred ages (due to their often large 
random and syst ematic uncertainties; see, e.g., 
Hillenbrand et all 2007). 

YSOs in their early phases are characterized 
by IR excesses due to the presence of a circum- 
stellar envelope and/or a disk that contains warm 
dust. The first IR classification in Class I, II and 
III sources, which is based on the IR spectral in- 



dex, was proposed bv lLadal (|l987t ). The discovery 
in the sub-mm of cold objects, younger than the 
IR sourc es, leads to the intr oduction of the Class 
sources (|Andre et al.lll993l ). Based on their mid- 
IR proper ties, YSOs are u sually classified as fol- 
lows (e.g. lAllen et al.|[2004 : Me geath et all 12004 ; 
Lada et al.l 120061 : I Allen eTall i20l)7h : Class ob- 



jects are protostars surrounded by the initial in- 
falling envelope, typically detected in the sub- 
millimeter and far-IR; Class I objects are evolved 
protostars surrounded by an envelope and a cir- 
cumstellar disk; Class II stars are objects with a 
circumstellar disk but lacking significant evidence 
for a dense envelope and finally Class III stars are 
objects that have lost their initial envelope and 
disk and show optical and near-IR colors consis- 
tent with normal photospheres. By providing very 
sensitive 3.6-70 /im data with good angular reso- 
lution, Spitzer observations have provided a much 
better census of the stellar populations in a num- 
ber of star-forming regions. However, Spitzer data 
have also complicated the YSO classification sys- 
tem, allowing detection of some Class objects at 
wavelengths as blue as 3.6 ^m. 

Our intent is to compare the X-ray properties of 
protostars - where to date there are relatively few 
studies - with the X-ray properties of the more 
evolved Class II and HI stars, for which the X- 
ray emission is well studied. Such a comparison is 
fundamental in order to understand the onset of 
X-ray emission in YSOs and its evolution in the 
early phases of star formation. 

Previous studie s of X-ray emission from Class 
O-I objects include lOzawa et al.l ([2005|), who stud- 
ied core F of the pOphiuchi cloud and found a 
high X-ray detection rate (64%) of Class I sources. 
They suggest that an evolutionary effect in tem- 
perature and extinction is present from Class I to 
Class HI stars, these latter having lower coronal 
temperatures and extinction. 

A s i milar result has also be en found bv lCiardino et al" 
(l2007l) . iGetman et all (|2007l) studied the IC 1396N 
cometary globule where they found 6 X-ray emit- 
ting protostars, including the Class O-I IRAS 
21391-1-5802 which they claim to be one of the 
youngest objects ever detected in X-rays. The 
comparison of X-ray emission from different evo- 
lutionary classes is, however, hampered by the 
difficulty of detecting X-rays from Class proto- 
stars. For example, no X-ray emission has been 
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found from Class star s in the Serpens, one of the 
most sensitive surveys ( Giardino et al.l[2007 ). even 
if important detections have been found in other 
surveys fe.g. iTsuboi et al.ll200ll: iHamaguchi et al 



20051: iTsuiimoto et al.ll2005l : iGetman et al.ll20071 ). 

The ONC is a prime target for this kind of in- 
vestigation. It is the nearest concentrated site of 
star formation with a very large number of known 
members covering all evolutionary stages. It is 
the best studied region at all wavelengths, in par- 
ticular, in the X-ray band. Many interesting re- 
sults have been obtained with the extraordinarily 
long Chandra observation performed by the COUP 
project. To understand the evolution of X-ray ac- 
tivity from the very first stages of the star forma- 
tion process, it is necessary to distinguish Class O-I 
objects from Class II stars. Such a classification, 
based on the shape of the spectral energy distri- 
bution (SED) is not available in the literature and 
requires complete and deep near and mid IR data. 

We present here the results of the X-ray study 
of YSOs in th e ONC classified usi ng the Spitzer 
IRAC catalog (jMegeath et al.ll2005l ). groundbased 
JHK data (Robberto et al. 2007, in preparation) 
and the deep photometric catalog obtained within 
the H ST Treasury Program (jRobberto et al 



2005a). Using the COUP data, we have stud- 
ied the X-ray properties of YSOs from their ini- 
tial protostellar phase (Class O-I objects) to the 
more evolved pre-main sequence (PMS) stage of 
no longer accreting and disk free stars (Class III 
objects). We have classified the ONC members 
using conservative criteria in order to minimize 
the contamination of the resulting samples. Be- 
cause of the large area coverage and the richness of 
the ONC, we have obtained an homogeneous and 
statistically significant sample of X-ray detected 
Class O-I objects. We are thus able, for the first 
time, to compute the X-ray luminosity function 
(XLF) of Class O-I stars and to compare it with 
the XLF of Class II and III stars. 

2. OBSERVATIONS 

We have assembled a deep photometric catalog 
of the ONC from the UV to 8 /xm obtained us- 
ing new optical/NIR observatio ns taken within 
the HST Treasury Program (jRobberto et al 



Fazio et al.ll2004) imaging. The optical photom- 
etry was derived using images taken with the 
Wide-Field Planetary C amera 2 (WFPC2) of HST 
(jRobberto et al.ll2005bl ) and with the Wide Field 
Imager (WFI) at the ESO 2.2 m Telescope (Rob- 
berto et al. 2007, in preparation). JHK pho- 
tometry was obtained from images taken with the 
Infrared Side Port Imager (ISPI) camera at the 
CTIO Blanco 4m telescope (Robberto et al. 2007, 
in preparation) . Spitzer IRAC photometry in the 
four standard filters ce ntered on 3.6, 4 5, 5.8 and 
8.0 /xm was obtained bv lMegeath et all (|2005l ). A 
detailed description of the optical-IR observations 
is given in Table[T] We note that the magnitude 
limits given in this table are typical values for 
each survey. However, due to the complex struc- 
ture of the nebulosity in the ONC (and due to 
varying integration times for different parts of the 
region), the sensitivity is not uniform. For ex- 
ample the limiting magnitude at 3.6 /im is about 
16 within most of the COUP FOV, while it is 
about 10 within a region of 1' in radius around 
the brightest stars of the Trapezium cluster. 

2.1. X-ray data 

X-ray data used in this work were obtained 
with the -838 ks COUP observation of the ONC, 
taken with the ACIS camera that covers a FOV 
of 17'xl7'. Details of data reduction and the 
derived catalog of 16 16 sources can be found in 
Getman et all j2005b[) . 



2005bl ). using a variety of telescopes and de- 
tectors, and Spitzer IRAC () Werner et al. 120041 : 



2.2. Optical data 

New optical photometric data t aken with the 
WFPC2 camera of HST were used bv lRobberto et al 
(|2005bl ) to derive the optical photometric catalog 
used in this work. This catalog contains 1754 se- 
lected point sources with at least a detection in 
the F814 band. 

The recent HST survey has been complemented 
with UBVI images taken with the WFI camera of 
the ESO 2.2 m Telescope. Details on data reduc- 
tion and photometric catalog derivation will be de- 
scribed in Robberto et al. (2007, in preparation). 
This catalog contains 2744 selected point sources 
with at least a detection in the I band. 

Additional optical data were ret rieved from 
l iterat ure including the catalogs of Hillenbrand 
()1997L hereafter H97), iRebull et all ()2000 h and 
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Table 1: Log of the new optical-IR observations used in this work. 



Nights 


Telescope 


Camera 


Filters 


Exp. time^ 


Mag. Lim.'^ 


Field of 


Reference 










[sec] 




view 




2006 May 


HST 


WFPC2 


F336W'= 


2x400 


F336W~20.0 


30' X 30' 


1 


2006 May 


HST 


WFPC2 


F439W'= 


80 


F439~20.7 


30' X 30' 


1 


2006 May 


HST 


WFPC2 


F656N'= 


400 


F656N~18.0 


30' X 30' 


1 


2006 May 


HST 


WFPC2 


F8UW 


10 


F814~20.0 


30' X 30' 


1 


2005 Jan. 1-2 


ES0 2.2m 


WFI 


U, B, V, Ic 


3, 30, 280 


Ic =^21.7 


35' X 34' 


2 


2005 Jan. 1 


CTIO Blanco 4m 


ISPI 


J, H, K, 


3, 30 


K, -18.7 


6x(10;25xl0;25) 


2 


2005 Jan. 2 


CTIO Blanco 4m 


ISPI 


J, H, K, 


3, 30 


K, -18.7 


8x (10^25x10^25) 


2 




Spitzer^ 


IRAC 


(3.6, 4.5, 5.8, 8.0) Atm 




[3.6]~16.8 


5.6 sq. degrees 


3 



"Exposure time per filter 

''Magnitude limit computed for stars with photometric errors smaller than 0.1 

=WFPC2 filters F336W, F439W, F656N, F814W corresponds to the standard filters U, B, Ha and I, respectively 

''Data part of the Spitzer/GTO Programs 43 and 50 taken in high dynamic range mode using 0.6 and 12 sec at each position. 



References. — fll lRobberto et al.l ll2005bl) : (2) Robberto et al. (2007) in preparation; (3') iMegeath et al.l ll2005t) 



Herbst et alj (|2002f) . 



HA[dc, 



Fig. 1. — The field of view (FOV) of the catalogs 
considered in this work are overplotted on the spa- 
tial distribution of the stars detected with IRAC. 
Note that the FOV of the IRAC (dashed line), 
RebuU et all (|2000l ) and 2MASS data (not indi- 
cated here) are larger than the field shown in this 
figure. 



2.3. JHK data 

The JHK photometric catalog has been ob- 
tained using images taken with the ISPI camera at 
the CTIO Blanco 4 m Telescope. Details on data 
reduction and photometry will be given in Rob- 
berto et al. (2007, in preparation). The survey 
includes a mosaic of 11 fields, each of 10!25x 10'25, 
centered around the Trapezium cluster. All fields 
were observed once using a 5 point dithering pat- 
tern. From the source catalog we excluded satu- 
rated objects, that is, those with a peak fiux in 
each band higher than 37000 counts and entries 
with errors in magnitudes larger than 0.2. After 
this selection, we were left with 17824 entries in 
the table. 

Finally, we constructed a consolidated catalog 
of individual objects by merging entries with po- 
sitions within 0"6, the radius at which spurious 
identifications start to dominate. For each object, 
we computed the mean coordinates and magni- 
tudes. Our final JHK photometric catalog con- 
tains 13 987 objects. 

Additional JHK data were retrieved from the 



2MASS catalog (|Cutri et al.ll2003l ). 



2.4. Spitzer IRAC data 

A description of the data, the reduction pro- 
cedure and the derived photometric catalog ob- 
ta ined from the Spitze r IRAC images is given 



Megeath et al 



( 20051) . The survey covers the 
Orion A and B molecular clouds for a total of 5.6 
sq. degrees with a sensitivity of 17.0, 16.0, 14.5 
and 13.5 mag in the 3.6 /im, 4.5 /im, 5.8 /im and 
8.0 /im bands, respectively. The uncertainties in 
the calibration are about 5%, but for a few stars 
are significantly worse due to the bright nebulos- 
ity in the 5.8 and 8.0 /im bands. Within the re- 
gion studied in this work, i.e. the COUP FOV, 
the catalog includes 1244 sources detected both at 
3.6 /im, 4.5 /im and 334 sources detected both at 
5.8 //m, 8.0 yum. The Spitzer Orion Survey also in- 
cluded MIPS observations at 24 //m that we do not 
use in the following analysis since they are satu- 
rated in the COUP region due to the high back- 
ground from the HII region. 
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2.5. Cross identifications of sources 

All the photometric catalogs obtained from 
the observations described in Table[T] were as- 
sembled in a single database. Common objects 
between two catalogs were found using the near- 
est neighbour method by first considering two 
sources as the same object if their positions agree 
to within a specified amount - l"-3"- depending 
on the absolute astrometric accuracy of the cata- 
logs. Matched sources were then used to compute 
the median and the standard deviation of the off- 
sets. These median values were used to remove the 
systematic offsets between each pair of matched 
catalogs, and new matches were performed, for all 
pairs of catalogs, using new tolerances computed 
using the standard deviation of the offsets. 

Using the same matching procedure we have 
cross-correlated these data w ith the literat ure cat- 
alogs ofiHi Uenbrand (,1997| ). iRebull et al.) (12000), 



2MASS fcutri et al| l20ol. iHerbst et al. [ (|2002), 
Muench et al.l ( 20021) a nd the X-ray source hst 



from the COUP Survey (jCetman et al.ll2005bf ) 

Finally we collected all data from all the cross- 
matched catalogs, including all objects detected in 
at least one catalog. 

The field of view (FOV) covered by these sur- 
veys are overplotted in Fig. [1] on the spatial distri- 
bution of the stars detected with IRAC around 
the ONC. Note th at the FOV of the IRAC, 
iRebull et all (|2000[ ) and 2MASS data (not in- 
dicated) are larger than the field shown in the 
figure. 

Since this work is focused on the X-ray prop- 
erties of the young stars in the ONC, only stars 
within the COUP FOV will be considered in the 
following analysis. In addition, photometric data 
with errors smaller 0.2 mag will be considered. 

3. YSO CLASSIFICATION 

In this Section we describe how we select ONC 
stars belonging to the three different evolutionary 
classes O-I, II and III. Our approach is to select 
conservative, even if incomplete, samples in order 
to be able to compare the X-ray properties of the 
stars in different evolutionary stages using samples 
that are as "uncontaminated" as possible. 



3.1. Class O-I protostars 

According to lLada fc Wilki"iie ( 19841 ) and Adams et al 



(|1987[ ). the first stage of YSO evolution is the in- 
fall phase around a central protostar undergoing 
accretion, characterized by strong emission in the 
sub-millimeter and far-infrared and weak emis- 
sion shortward of 24 ^m. These objects, indicated 
as Class 0, can have outflows and are defined as 
protostars with half or more of their mass in the 
envelope. More evolved protostars are the Class I 
objects with both a disk and envelope surrounding 
the central star. Both these classes of protostars 
are characterized by a steeply rising SED in the 
infrared range, although pure Class objects are 
expected to be fa inter at near and mid-IR bands 



(| Allen et al.l 120071 ) 



Since our data extend only to 8 /im, we cannot 
distinguish between Class and Class I stars and 
we will combine them into one group as Class 0- 
I protostars; we classify them using the following 
criterion: we computed log(A-F'A) over the wave- 
lengths [Ai,...A5] = [2.2, 3.6, 4.5, 5.8, 8.0] Aim, i.e. 
the K-band and the four IRAC bands with the 
Lada et al. I (|2006l) flux conversion factors that are 



in agreement with the Reach et al. I (|2005l) values; 
for each star we then computed the slopes for each 
pair of adjacent bands 

log(A,+iFA,+ J - log(A,f^A; 



for 



1,4; 



logAi+i - logAi 

we consider Class O-I objects those for which the 
slopes are all larger than 0, i.e. objects with ris- 
ing SEDs. This criterion is more conservative 
than using the single slope (usually indicated as 
spectral index) computed from the linear fit of 
the log(Ai^A) function over the whole wavelength 
range, because it only selects monotonically in- 
creasing SEDs. Moreover, it is more robust than 
those based on only a color-color diagram (CCD) 
since it is based simultaneously on all the known 
magnitudes. 

Our sample of Class O-I objects includes all the 
objects for which we have at least three of the 
five magnitudes and therefore at least two slopes. 
This allows us to include objects detected only in 
the K, 3.6 fiui and 4.5 /im bands but not at 5.8 & 
8.0 /im; because of the bright nebular background 
and the lower sensitivity at 8 /im, many objects 
are not detected in these two bands (roughly 44% 
of the Class O-I sample). 
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Using this criterion, we end up with 41 poten- 
tial Cla ss O-I obje cts in the COUP field of view 
ICetm a n et alJ [2005b) . Eighteen of them have 
a detection in at least one of the optical UBVI 
bands. Since protostars are expected to be heavily 
embedded and therefore highly reddened objects, 
we do not expect to detect them in the visible. Of 
these 18 peculiar objects, 5 (COUP # 693, 747, 
826, 948, 1011) a re in the list of t he X- ray detected 
proplyds given in lKastner et al. I (|2005h . From a vi- 
sual inspection of ACS/HST images, we conclude 
that all 18 of these objects have a counterpart 
in the HST optical bandpasses, while the other 
23 objects are not detected at these bands. We 
eliminate these 18 objects from our list of candi- 
date Class O-I stars and we consider the remaining 
23 stars as Class 0-Ia objects to distinguish them 
from the additional sample of candidate Class 0- 
I described below. The SEDs of the 23 selected 
Class 0- la stars are plotted in Fig. [51 The pho- 
tometry is given in Table[5] where Column 1 gives 
the sequential number in our catalog. Column 2 is 
the COUP identification number. Columns 3 and 
4 are the J2000 coordinates. Columns 5 to 11 give 
the JHK and IRAC magnitudes. Of these 23 ob- 
jects, 10 (~40%) have an X-ray counterparts in the 
COUP data. One of t hese (COUP # 702) is in the 



Getman et all (|2005ar ) hst of nonflaring COUP X- 



ray sources without optical and NIR counterpart 
and was flagged by them as a probable extragalac- 
tic object (EG). 

Fig. [3] shows the color-color diagrams (CCDs) 
obtained using the K and IRAC magnitudes of 
all the stars in the COUP field of view. The 
plot on the left panel includes more stars than 
those on the central and right ones because it does 
not require detections at 5.8 and 8.0 /xm. The 
limited sensitivity of the 8.0 /im channel means 
that few objects with zero color (consistent with 
photospheric emission) appear in the right panel. 
The selected Class 0-Ia objects with increasing 
SED, indicated in Fig. [3] by squares, show [3.6]- 
[4.5]>0.7 and [5.8]-[8.0]>l.l. The objects with 
[3.6]-[4.5]> 0.7 and [5.81-[8.0]<1. 1 that are usually 
also classifi ed a s Cla ss O-I stars ( Allen et al.ll200^ 
iMe geath et 'aLl [200l . show a decreasing or flat 
SED at wavelengths longer than 4.5 /im. The spec- 
tral energy distributions of protosta rs in this spec- 



reduce the extinction toward the poles of the 
protostellar envelopes allowing radiation to leak 
and scatter out. Consequently, protostars can 
also exhibit decreasing slopes, particularly be- 
tween 4.5 and 8.0 fun , depending on their inclina- 
tion ( Whitnev et al.l 2003). Furthermore, highly 
reddened Class II objects can have increasing 
slopes; however, these slopes will generally flat- 
ten or decline for A > 4.5/xm due to the flattening 
of the redd ening law and the bro ad silicate feature 
at 9.7 /xm (jPlahertv et al.ll2007t ). Thus, although 
the IRAC data can enhance the identification of 
protostars by extending the SEDs to 8 microns; es- 
tablishing a sample of protostars uncontaminated 
by reddened Class II objects is still difficult with- 
out photometry at wavelengths A > 8 /im. 

To take into account these cases, we also include 
the objects with a decreasing SED for A > 4.5/im 
in our Class O-I sample. There are 28 such ob- 
jects in total. After rejecting the six objects with 
optical counterparts in the HST imaging, we are 
left with 22 objects, which we refer to hereafter 
as Class 0- lb sources. The SEDs of the Class O-Ib 
objects are shown in Fig.lH while the photome- 
try is given in Tabled analog to TableO Among 
these objects, 18 (-80%) are in the COUP hst 
of X-ray detection, including COUP source 1565, 
which was classified a s an extragalactic object by 



Getman et al.l gOOM) 



If we consider both Class 0-Ia and Class 0-Ib 
samples, our list of X-ray detected Class O-I ob- 
jects (10-1-18) is one of the largest sample of 
Class O-I stars with evidence of X-ray e mission 
(jCetman et al.ll2007t Iciardino et al.ll2007^ . 

Our selected Class O-I objects fall i n the region 



of the enve lope models in d icated in lAllen et al 
(|2004l ) and iMegeath et"all (|2004f ). There are a 



tral r egime are complex as r otation (IKenvon et al 



19931 ) and outflow cavities (jWhitnev et al.l l2003f ) 



small number of objects falling in the Class O-I re- 
gions of the CCD which were not selected with 
our criterion because either they have an optical 
counterpart (18-1-6 in total) or their SED band- 
to-band slopes do not increase as required by our 
criterion. Whereas some of these rejected objects 
could be Class O-I stars, by adopting a conservative 
selection criterion, we have minimized the adverse 
effect of sample contamination for our comparison 
of X-ray properties. 

A careful examination of the near-IR and IRAC 
colors of our Class O-I candidates (see Figs.|3| and 
[5]) suggests to us that the Class O-Ia's (on average) 
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Fig. 2. — SED of the 23 Class 0-Ia candidate objects with a rising SED from K to S.Ofxm. The sequential 
number in our catalog and the spectral index computed using the available IRAC magnitudes are also 
indicated. 
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Fig. 3. — CCD of the stars in the COUP field of view (dots) created using the K and the Spitzer IRAC 
magnitudes. The total number of objects in each diagram is reported in the three panels. Squares and 
large filled circles represent the Class 0-Ia and 0-Ib objects, respectively, while X symbols indicate those with 
an X-ray counterpart; asterisks indicate objects rejected becau se they have an optic al counterpart. The 
reddening vectors computed using the IR reddening law given in iFlahertv et al.l (|200l1 ) for Orion A are also 
plotted. 



are inconsistent with their being heavily reddened 
Class II's. In particular, the bottom-right panel of 
Figure 5 shows that, using the Orion A reddening 
law of Flaherty, the Class 0-Ia's are too red in [4.5]- 
[5.8] to be heavily reddened Class IPs - and in any 
case would require reddenings of Ay > 60. The 
bottom left panel of Fig.lH] also indicates that the 
Class 0-Ia's, on average, would need very large red- 
denings (Ay ~ 50) for them to be Class II sources. 
Such heavily embedded objects would likely be 
very young in any case - and we therefore believe 
it is more likely that they are, indeed. Class or 
Class I objects. By contrast, the Class 0- lb objects 
are located in all of the diagrams in Figs.[3]and[5] 
in positions consistent with being moderately red- 
dened (Ay ~ 30 on average) Class IPs. We suspect 
the Class 0-Ib's are a mixture of true ClassO or I 
objects, and heavily reddened Class II's - without 
more data we cannot clean the sample any further 
than we have. 

Another possible worry is that our Class O-I 
candidate sample includes systems seen edge-on 
in or around the HII region. In these cases, rela- 
tively evolved Class II objects may display a SED 
peaking a t mid-IR and far -IR wavelengths, as sug- 
gested by Robberto et al.l (2002). However, since 
edge-on disks should be randomly distributed in 



the nebula, clustering in the dark regions of the se- 
lected objects (see Fig. [5^) strongly suggests that 
they are real protostars. 

3.2. Class II stars 



Following iHartmann et ah ( 2005 ). Class II ob- 
jects within the COUP field of view were selected 
using the condition 0.2 < [3.6] - [4.5] < 0.7 and 
0.6 < [5.8] - [8.0] < 1.1. These criteria are 
slightly more r estric tive than those adopted by 
Megeath et al. ( 20041 ) in order to attempt to create 
as "pure" a set of Class IPs as possible. With our 
selection criteria, we may lose reddened Class II 
objects, but assuming that X-ray properties of 
Class II stars are independent of their reddening, 
our sample should be representative of the Class II 
population. 

Using this criterion, we find 148 Class II stars 
with NIR photometry given in Table[71 analogous 
to TableEl 130 of these objects (about 88%) are 
X-ray detected. One of them, the CO UP source 
1401, is listed bv ICetman et aL I (|2005al) as a flar- 
ing source and classified as an Orion Molecular 
Core (OMC) member, whereas another Class II 
object, COUP source 751, was classified by the 
same work as extragalactic object. Among the 130 
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Fig. 4. — SED of the 22 Class 0-Ib candidate objects with a rising SED from K to 4.5 /xm. The sequential 
number in our catalog and the spectral index computed using the available IRAC magnitudes are also 
indicated. 
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X-r ay detected Class I I stars with a counterpart in 
the Hillenbrandl (1997) catalog, 83 have a measure 
of the equivalent width of the Call 8542A line 
W(Call). Among these, 28 have VF(CaI I ) < -1, 



which is the limit used bv lFlaccomio et alj (|2003al ) 
to classify strong accretors. Therefore we have 
about 55 objects (i.e. ~66%) that have a near-IR 
bright inner disk but relatively low accretion rate. 
Of the 148 Class II stars we have identified, we 
will consider in the following analysis only the 87 
X-ray detected and the 9 X-ra y undetected stars 
for which masses ar e known ( Hillenbrandl Il997 ; 
Getman et al.l[2005lJ l. We retain only these stars 
in order to study the X-ray properties as a func- 
tion of stellar mass. 

3.3. Class III stars 

Class III objects are PMS stars that do not 
show IR excesses and therefore were selected as 
those having IRAC colors near zero. In order to 
include all the stars not detected in the longer 
wavelength IRAC bands, we selected all the stars 
with X — [3.6] < 0.5 and rejected those with [3.6] — 
[4.5] > 0.2 or [4.5] - [5.8] > 0.2 or [5.8] - [8.0] > 0.2, 
when these latter colors were available, i.e. objects 
with IR excesses. Stars undetected in the IRAC 
bands at wavelengths longer than 3.6 ^m were in- 
cluded since their non-detection in these bands is 
compatible with purely photospheric SEDs. 

Using these criteria, we selected 205 objects, 
with photometry given in TablelU 160 of them 
(about 80%) have an X-ray counterp art. Amon; 
these latter, 150 are in the list of iHillenbran 
()l997l) . In particular, 130 of them have a proper 
motion membership probability larger than 90%, 
6 have a membership probability equal to and 
the other 14 have no membership probabilities. 
For the remainin g 10 stars w i thout a counter- 
part in the list of Hillenbrandl (1997), X-ray de- 
tection is the only known membership criterion. 
We note that none of these 10 objects is either in 
th e list of flaring or in t hat of non-flaring sources 



of Getman et al 



(|2005bl ). Therefore the member- 
ship class for these objects remains uncertain. 

Of the 45 objects without an X-ray detection, 
18 have a counterpart in the Hillenbran d (199_^ 
catalog but only 5 of these have a proper motion 
membership probability larger than 90%. 

The remaining 40 objects without a mem- 



bership criterion could be either Class III stars 
without X-ray activity, or field stars. Of the 
205 Class HI stars initially selected, we will con- 
sider the 103 X- ray detected stars for which the 
mass is known (Getman et al. I l2005bh and with 
Hillenbrand ( 19971 ) membership probability larger 
than 90%; in addition, from the sample of X- 
ray undetected Class HI candidates, we will retain 
only the 2 objects with membership probability 
larger tha n 90% and for whi ch a mass has been 
derived bv lHillenbrai^ (|l997l ). 

Fig. ID shows the selected stars in the IRAC 
CCD and color-magnitude diagrams. These di- 
agrams show that our samples of YSOs most 
likely do not include AGN or star-forming galax- 
ies; that is, they do not include objects falling in 
the wedge-shaped regions outlined in the upper 
panels of Fig.[5l Star-forming galaxies are usu- 
ally dominated by PAH features in the 5.8 and 
8.0 /im and therefore have [5.8] - [8.0] > 1.0 but 
[3.6] — [4.5] < 0.3, that corresponds to the regions 
indicated by solid fines in the t wo upper panels 
of Fig.lSl (|Gutermuth et"al]l2007b . AGN, instead, 
have [3.6] — [4.5] and [5.8] — [8.0] colors very sim- 
ilar to the Class II YSOs, but are typically much 
fainter than YSOs, having, for example, [4.5] mag- 
nitudes (Megeath et al.l 120051 ). much fainter than 
those of the selected YSOs (see bottom-left panel 
in Fig.Ej). 

4. SPATIAL DISTRIBUTION 

The spatial distribution of YSOs within the 
molecular cloud from which they have formed may 
give important clues for understanding the star 
formation process. Using the classification defined 
in the previous section, we have analyzed the spa- 
tial distribution of the YSOs in the ONC, as shown 
in Fig. [6^ where we plot the position of all 2590 
objects in our catalog within the COUP field of 
view that have been detected in at least two of 
the UBVIJHK and IRAC bands and the Class O-I 
candidates. Stars of Class II and HI arc plotted in 
Fig.[6l3 and c, respectively. 

Figure!?] shows t he velocity i ntegr ated ^'^CO 



map obtained by iBallv et al.l (|l987l ) and the 



Spitzer IRAC image at 3.6 /im where the dif- 
ferent classes selected in this work have been 
overplotted. The upper green box indicates a 
40"x50" region around the Becklin-Neugebauer 
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Fig. 5. — IRAC CCD and color-magnitude diagrams of the different populations in the ONC. Dots indicate 
all the stars in the COUP field of view with detections at the bands used in each panel. Squares and large 
filled circles are the Class 0-Ia and 0-Ib objects, respectively, plus symbols are the stars classified as Class II 
and diamonds are th e Class III stars. The reddening vector corresponding to Ay = 30, computed using the 
Flahertv et aD (|2007t ) reddening law, is plotted on each panel. The solid lines in the above panels separate 



the region dominated by YSO's from the region expected to be dominated by galaxies. 
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Fig. 6. — Spatial distribution of all the objects in our catalog within the COUP field of view that have been 
detected in at least two of the UBVIJHK and IRAC bands. Squares and large filled circles indicate Class 0-Ia 
and 0-Ib objects, plus symbols the Class II stars and diamonds the Class III objects. In all the panels, X-ray 
detected objects are indicated with an X symbol. 



object and the Kleinmann nebula, also indicated 
as BN-KL (RA2ooo=05''35™14':16, Dec2ooo=- 
05°22'21.5"). The bottom green box indicates a 
60" X 75" region aroimd the Orion Molecular Cloud 
1 (OMC-IS) (RA2ooo=05''35™14^50, Dec2ooo=- 
05°25^49.0"). Th ese boxes show the limits of the 
areas studied in lCrosso et al. ( 2005h . 

Most of the Class O-I stars are concentrated in 
the central and north-east parts of the COUP field 
of view along a region running north-south from 
the Trapezium that follows the north-south dis- 
tribution of the ^^CO emission, i.e. of the dense 
star-forming gas. The distribution of our selected 
Class O-I objects is biased by our NIR observations 
which are more sensitive towards the north of the 
molecular ridge where the nebular contribution is 
lower. The lack of Class O-I objects in the BN-KL 
and OMC-ls regions is due to the bright nebulos- 
ity, as is evident from the image at 3.6 /im, that 
causes a decrease of sensitivity in the central re- 
gion. For this reason, the spatial distribution of 
our Class O-I objects cannot be directly comp ared 
with that obtained bv lLada et aD (|200(ll2004[ l. al- 
though in fact the two spatial distributions are 
quite similar. 

Figure s[6land[7lshow th e same structure already 
found in Feigelson et al.l |2005) using the hghtly 
X-ray absorbed stars: they interpreted this distri- 
bution as due to a process of violent relaxation. 



Note that, the latter study examined the spatial 
distributions stratified by absorption, while here 
the stars are stratified by individual YSO evolu- 
tionary classification. The features of the spatial 
distribution visible in Figs. [6] and [7] include: the 
OMC-IS core, BN/KL core and northern exten- 
sion which are part of the larger filament along the 
middle of the Orion A clo ud and finally th e East - 
West asymmetry that in iFeigelson et al. (|2005l) 
was attributed to non-equilibrium dynamics. The 
latter only shows up in Class III systems. 

The spatial distributions of Class II and III 
stars are also in ag reement with the s patial dis- 
tri bution found by Hillenbrand! ( 1997 ) who, us- 
ing D'Antona fc Mazzitellil ( 1994 ) isochrone ages, 
found that youngest stars are clustered toward the 
cluster center while slightly older stars are more 
widely distributed. The elongated asymmetric dis- 
tribution of our candidate Class O-I stars and the 
large-scale EW asymme try are consistent with the 
distrib ution derived by iHillenbrand fc Hartmanii 
(jl998[ ) who concluded that the structure of the 
ONC is more the result of the different environ- 
mental conditions rather than the result of dynam- 
ical evolution. 

There is a lower percentage of Class II objects in 
the central region, but this could be due to the re- 
liance on the 8.0 /xm detections for the selection of 
Class II stars. Our spatial distribution are, there- 
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Fig. 7. — (Left) Velocity integrated ^'^CO emissivity from iBallv et"al (1987) in the region around the 
Trapezium cluster and (right) Spitzer IRAC image at 3.6 /xm in the same region. Both images are shown 
with inverted color map. Yellow squares and light blue circles indicate Class 0- la and 0-Ib objects, green plus 
symbols are Class II stars while red diamonds are Class III stars. The 40"x50"upper green box indicates the 
BN-KL region while the 60"x75"bottom green box indicates the Orion Molecular Cloud 1 South (OMC-IS). 



fore, strongly dependent on the sensitivity of our 
NIR data which is quite limited toward the central 
region, where the contamination by the bright OB 
stars and nebular emission is very high. 

5. EVOLUTION OF X-RAY ACTIVITY 

The very deep X-ray COUP data together with 
the relatively large sample of Class O-I objects se- 
lected with the Spitzer observations allow us, for 
the first time, to analyze and compare, in a statis- 
tically significant way, the X-ray properties of the 
YSOs in their initial protostellar phase with those 
of more evolved PMS stars. 

As already mentioned in Section[31 we find that 
the fraction of Class 0-Ia objects detected in the X- 
rays is almost 40% (10/23), that of Class 0-Ib and 
Class III is about 80% (18/22 and 160/205, respec- 
tively) while that of Class II objects is about 88% ( 
130/148). The analysis described in this section 



aims at understanding whether the X-ray emis- 
sion can be related to the IR properties of YSOs, 
which we take as signposts for the different phases 
of protostellar and PMS evolution. 

5.1. X-ray luminosities 

We use the results of the spec tral analysis pre- 
sented by ICetman et al. (|2005d) to compare the 
X-ray properties of the different populations of 
YSOs in the ONC. For the X-ray luminosities of 
detected stars we adopt the absorption corrected 
luminosities in the total 0.5-8 .0 keV band, Lt^c, 

( 2005bl) . 



from Table 8 of Getman et al 



In order to compare unbiased luminosity func- 
tions of the different classes, we derive upper lim- 
its to the X-ray luminosities of the X-ray un- 
detected objects. Upper limits to the photon 
count rates of the 13-1-4 undetected Class 0-Ia and 
Class 0-Ib objects were calculated with PWDe- 
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tect ( Damiani et al.|[l997al lbl) consistently with the 
source detection procedure o f the other COUP 
sources ( Getman et al. adopting a detec- 

tion threshold significance of 5 a. The upper limits 
for the 9 undetected Class II stars and for the 2 un- 
detected Class III stars were taken from Table 11 



of iGetman et alJ (|2005bl ). 

In order to convert upper limits to the count 
rates into unabsorbed X-ray luminosities, we de- 
rive an extinction dependent conversion factor. 
The CCDs in Fig. [3] indicate that X-ray detected 
and undetected Class O-I objects have similar red- 
dening. We therefore computed the conversion 
factors of the X-ray detected Class O-I objects as 
the ratio between the unabsorbed X-ray luminosi- 
ties (Lt,c) and the count rates taking into account 
the PSF fraction within the extracted area and 
the effective exposure time at the source position 
( Getman et al. lioOSbh . and we plotted these val- 
ues as a function of the NjJll values. We assume for 
the X-ray undetected Class O-I objects a conver- 
sion factor computed as the median value of the 
conversion factors of the X-ray detected Class 0- 
I objects, which corresponds for these stars to 
A^H — 4.5 X lO^^cm"^. However, in order to check 
whether our results are robust with respect to the 
assumed Njj of undetected objects, we also consid- 
ered the conversion factors corresponding to the 
highest and minimum conversion factors of the X- 
ray detected Class O-I objectsU We used the same 
procedure to estimate the upper limits to the X- 
ray luminosities in the 2.0-8.0 keV hard band. 

The X-ray properties of the Class O- I COUP 
X-ray sources, taken from lCetman et aL ( 2005bl ). 
are given in Table[9l We provide sequential num- 
bers from our catalog, sequential COUP numbers, 
background-corrected counts in the 0.5-8.0 keV 
band, the fractions of the PSF within the source 
extraction areas, effective exposure times, the hy- 
drogen column densities obtained from the spec- 
tral fit, plasma temperatures, emission measures 
and finally the absorption-corrected X-ray lumi- 
nosities in the 2.0-8.0 keV hard band, Iju.c, and in 
the 0.5-8.0 keV total band, Lt,c- 



^ The H column densities were derived by I Getman et al.l 
l l2005bl ) from the spectral analysis of the COUP sources. 

^We used the median conversion factor to derive the X-ray 
luminosity (log Lt,c = 28.52 erg/s) o f the COUP source 
1197 for which iGetman et al.l l|2005bl ) did not perform a 
spectral fit because of the low counting statistics. 



Values computed for the X-ray undetected 
Class O-I stars are given in TablefTOl We list se- 
quential numbers from our catalog, the celestial 
coordinates from the Spitzer catalog, upper limits 
to the source counts computed with PWDetect 
in the hard and total bands, the effective Chan- 
dra exposure time at the object positions and, 
finally, the upper limits to the unabsorbed X-ray 
luminosities in the hard and total energy bands. 

X-ray properties of X-ray detected Class I I and 



( 2005bh . 



are 



III stars, taken from Getma n et al. 
given in Tables [Tl] and [12] where column names 
are analogous to those given in Tabled 

The X-ray luminosities of X-ray undetected 
Class II and III objects were derived, as for Class 0- 
I stars, using the median conversion factor of the 
X-ray detected stars in the same class. This as- 
sumption is based on the fact that the NIR colors 
of the X-ray undetected Class II stars indicate that 
these objects are not more absorbed than those 
detected with COUP. In addition, the conversion 
factors corresponding to the lowest and highest ob- 
served Nh values were also considered (10^^ and 
4.0 xlO^^ cm~2 for the Class II stars and 0.5x10^^ 
and 2.0 x 10^^ cm^^ for the Class III stars). Values 
computed for the X-ray undetected Class II and III 
stars are given in TablefTSl analogous to TablefTOl 

5.2. X-ray luminosity functions 

We now compare the X-ray luminosity func- 
tions of the stars of the four considered NIR 
classes. In order to take into account nondetec- 
tions, XLFs were derived using the Kaplan-Meier 
maximum-likelihood estimator. For Class II and 
III sources, we considered separately stars in the 
three mass ranges used by IPreibisch &: Feigelson 
()2005l ). Fig. [8] shows the comparison of the XLFs 
of the four populations. The red and green lines 
indicate the XLF of Class III and II objects, re- 
spectively, with masses in the given range, while 
the black and blue lines are the XLFs of Class 0- 
la and 0-Ib stars, for which masses are unknown 
and are therefore the same in the three panels. 
Lines of the same color in each panel show the 
XLF obtained using the three different count rate 
to luminosity conversion factors adopted for the 
X-ray non detections of each subsample and give 
an idea of the uncertainty of the XLFs. We re- 
moved from this analysis the four objects of our 
sample of Class 0-Ia protostars without detections 
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Fig. 8. — Maximum likelihood cumulative XLF in the total (0.2-8.0) keV energy band of the four NIR 
classes. The three panels refer to Class II and Class III stars in three different mass ranges. The red and 
green lines indicate the XLF of Class III and II objects, respectively, with masses in the given range, while 
the black lines are the XLF of Class O-I stars, for which masses are unknown and are therefore the same in 
the three panels. Lines of the same color in each panel show the XLF computed using the three different 
luminosity conversion factors for the X-ray non detections. Thick lines indicate the XLF used to perform 
the statistical tests (see text). Note that the distributions do not reach the f=100% because of the presence 
of upper limits. 




Fig. 9. — Maximum likelihood cumulative XLF in the hard (2.0-8.0) keV energy band of the four NIR classes. 
The figure is analogous to Fig. [8] 



at 4.5 /im. Since this band is one of the most sen- 
sitive bands and the least contaminated by nebu- 
losity, these objects could be PAH contaminated 
sources or even knots of nebulosity. By discarding 
these objects, we are more confident to consider 
a sample of protostar candidates as less contami- 
nated as possible. The figure shows that for the 



two highest mass ranges, the X-ray luminosities 
of Class III stars are highly elevated with respect 
to those of Class II stars and 0-Ib objects. These 
latter two classes have very similar XLFs and are 
significantly brighter than Class 0-Ia stars. 

Because the Class 0-Ia and Class 0- lb stars 
in our sample are heavily absorbed, a note of 
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caution is needed. Since absorption precludes 
the observation of the soft part of the intrin- 
sic source X-ray spectra, we cannot rule out 
the presence of soft, completely absorbed emis- 
sion from Class O-I stars. Nevertheless, we know 
from previous studies of CTTS stars that a 
very soft component, kT<0.3keV, in the few 
cases in whi ch it is observed, is mainl y due 
to ac c retion ([Kastner et al. J 20021: S chmit t et alJ 



20051: iGiidel et all l2007t lArgiroffi et all \2m 



Huenemoerder et al. 2007: Grosso et al 



while harder components are due to coronal-like 
emission. This is in agreement with the previ- 
ous finding for Class II stars in other regions (e.g. 
Stelzer fc Schmittll2004l ). where the X-ray spectra 
of Class IPs were found to include an additional 
softer emission component compared to WTT's. 
Our Orion stars are too heavily absorbed to detect 
this soft component. 

Coronal emission, as analyzed with resolution 
spectra from Charge Coupled Devices, can usu- 
ally be decomposed into so fter ^0.8-1.0 keV an d 
harder ~2keV components ( Preibisch et al.ll2005[ ). 
The soft coronal component, like the accretion 
component, might remain unobserved in heavily 
extincted sources such as our Class O-I stars and 
thus the X-ray luminosities in the total band can 
underestimate the true emission in such objects. 
In addition, magnetic flare components can be 
strong around 0.5-2.0 keV, and dominate the X- 
ray flux. For these reasons, we also considered the 
comparison of the XLFs of different classes, by re- 
stricting the analysis to the X-ray luminosities in 
the hard (2.0-8.0) keV energy band. We note that 
we recomputed with PWDetect the upper limits 
in the hard band for a total of 31 undetected ob- 
jects, assuming a threshold of 4.8(t estimated from 
the background in this band. Among these ob- 
jects, 12 are flagged by the PWDetect code as af- 
fected by a X-ray detected source nearby, which 
therefore raises the upper limit count rate. In two 
cases out of 12, however, the asymmetry of the 
Chandra PSF at large off-axis angles is such as to 
make PWDetect clearly overestimate the contri- 
bution of the detected source to the computed up- 
per limit at the desired position, which is instead 
dominated by the background; in these instances 
we computed the upper limit from the background 
alone. The results of this analysis are shown in 
Fig. [HI analogous to Fig. [8] but for the hard energy 



band. 

To check if these results are statistically signif- 
icant, we performed statistical tests developed for 
"survival analysis" , suited for the analysis of cen- 
sored datasets. The probabilities that the X-ray 
luminosities of the stars in two samples are drawn 
from the same parent distribu tion were computed 
with the ASUKvd package (IPei gelson fc Nelson 
The results of this analysis are given in 



Tables[2]and[3]where, for each pair of samples, we 
list the probabilities computed using the five tests 
of the ASURV package. We note that the tests 
were performed in the most conservative way, i.e. 
by considering for each pair of samples the nearest 
XLFs among those computed with different con- 
version factors for the upper limits. The tests show 
that Class 0-Ia and 0-Ib stars have statistically dif- 
ferent XLFs both in the hard and total bands, at a 
significance level larger than 99%; Class 0-Ia and 

II are different at a significance level of 99% for 
the [0.9-1.2] and [0.5-0.9] M© mass ranges. The 
difference is marginal, both in the total and hard 
band, for the [0.1-0.5] M© mass range. The dif- 
ference between Class II and III objects is quite 
marginal (between 85% and 95%) for the highest 
mass range, while for lower mass stars. Class II and 

III stars do not show any significant difference. 
These results clearly indicate an evolution of 

the X-ray activity from Class 0-Ia, the least X-ray 
luminous objects, up to Class III, which show the 
highest X-ray luminosities, at least for stars with 
mass larger than 0.9 M©. Although the difference 
in the intrinsic X-ray activity between Class II and 
III stars is quite marginal, it is in agreement with 
previous results indicating that the X-ray lumi- 
nosities of accretors in th e ONC are smaller than 
those of non accret ors ( Flaccomio et al. l2003bl: 
Preibisch et al. Iliooi). Note, however, that these 



results cannot be compared directly because of the 
different criteria adopted to select the two popu- 
lations. 

The most significant and novel result is the dif- 
ference in X-ray activity level between Class 0-Ia 
and Class II and between Class 0-Ia and Class 0- 
Ib sources. Our results are even more robust if 
we include the four objects without detections at 
4.5 fim. 



^Astronomy Survival Analysis available from the StatCodes 
|http://www.astro.psu.edu/statcodes| 
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Table 2: Probabilities of correlation between the XLFs in the total band (0.2-8.0) keV of two samples, 
computed with the ASURV package (see Section [5T2l) . Pi and P2 are the probabilities computed using the 
Gehan generalized Wilcoxon test with permutation and hypergeometric variance, respectively, P3 is the 
probability from the Logrank test while P4 and P5 are the probabilities from the Peto and Peto and from 
the Peto and Prentice generalized Wilcoxon tests, respectively. 



Mass 


Pi 


P2 


P3 


P4 


P5 


Sample 1 


Sample 2 


[M 


3] 


[%] 


[%] 


[%] 


[%] 


[%] 










0.23 


0.21 


3.13 


0.36 


0.25 


Class 0-Ia 


Class 0-Ib 


[ 0.9- 


1.2 ] 


0.35 


0.07 


0.03 


0.46 


0.45 


Class 0-Ia 


Class II 


[ 0.9- 


1.2 ] 


94.75 


94.99 


62.55 


97.88 


97.73 


Class 0-Ib 


Class II 


[ 0.9- 


1.2 ] 


6.79 


5.10 


16.15 


16.15 




Class II 


Class III 


[ 0.5- 


0.9 ] 


0.04 


0.02 


0.06 


0.05 


0.03 


Class 0-Ia 


Class II 


[ 0.5- 


0.9 ] 


35.94 


35.60 


58.58 


36.77 


36.64 


Class 0-Ib 


Class II 


[ 0.5- 


0.9 ] 


33.13 


32.21 


14.63 


31.51 


32.68 


Class II 


Class III 


[ 0.1- 


0.5 ] 


5.28 


5.64 


23.77 


6.04 


5.53 


Class 0-Ia 


Class II 


[ 0.1- 


0.5 ] 


3.41 


2.02 


1.52 


3.43 


2.88 


Class 0-Ib 


Class II 


[ 0.1- 


0.5 ] 


86.08 


86.09 


40.72 


80.37 


80.34 


Class II 


Class III 



Table 3: Probabilities of correlation between the XLFs in the hard band (2.0-8.0) keV of two samples, 
computed with the ASURV package (see Section [5T2|) . The probabilities are those described in Table[2l 



Mass 


Pi 


P2 


P3 


P4 


P5 


Sample 1 


Sample 2 


[M 


3] 


[%] 


[%] 


[%] 


[%] 


[%] 










0.36 


0.36 


2.99 


0.53 


0.47 


Class 0-Ia 


Class 0-Ib 


[ 0.9- 


1.2 ] 


0.57 


0.15 


0.08 


0.72 


0.73 


Class 0-Ia 


Class II 


[ 0.9- 


1.2 ] 


94.84 


95.04 


46.15 


93.70 


94.39 


Class 0-Ib 


Class II 


[ 0.9- 


1.2 ] 


6.54 


4.77 


9.37 


9.37 




Class II 


Class III 


[ 0.5- 


0.9 ] 


0.08 


0.05 


0.89 


0.12 


0.08 


Class 0-Ia 


Class II 


[ 0.5- 


0.9 ] 


28.18 


27.57 


56.83 


29.37 


28.97 


Class 0-Ib 


Class II 


[ 0.5- 


0.9 ] 


62.80 


62.59 


34.13 


61.58 


62.47 


Class II 


Class III 


[ 0.1- 


0.5 ] 


14.23 


14.57 


45.40 


16.79 


16.58 


Class 0-Ia 


Class II 


[ 0.1- 


0.5 ] 


3.62 


2.19 


2.05 


3.70 


3.02 


Class 0-Ib 


Class II 


[ 0.1- 


0.5 ] 


46.63 


46.63 


91.93 


50.19 


50.23 


Class II 


Class III 
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Fig. 10. — Absorption corrected X-ray luminosity 
of Class 0-Ia stars (squares) and Class 0-Ib (circles) 
as a function of the [5.8]-[8.0] (panel a) and [3.6]- 
[4.5] (panel b) colors. Upper limits are indicated 
by arrows. 



Based on our selection criterion, Class 0- la ob- 
jects are protostars with SEDs increasing up to 
8.0 /im, while Class 0-Ib are protostars with SEDs 
increasing up to 4.5 /im and decreasing at longer 
wavelengths. Fig.fTUb shows the absorption cor- 
rected X-ray luminosities plotted vs. the [5.8]-[8.0] 
color, that discriminates Class 0-Ia stars (squares) 
from Class 0-Ib ones (filled circles), the former be- 
ing redder than [5.8]-[8.0]~l.l and the latter bluer 
than this value. This plot indicates a trend of 
increasing X-ray luminosity with decreasing [5.8]- 
[8.0] color. Fig.flOb shows the same luminosities 
as a function of the [3.6]-[4.5] color. While this lat- 
ter color does not discriminate between Class 0-Ia 
and 0-Ib objects, it clearly confirms that, from the 
X-ray point of view, Class 0-Ia and 0-Ib objects 
belong to two different populations. This again 
suggests us that the contamination by reddened 
Class II objects can be very high for Class 0-Ib ob- 
jects, while the Class 0-Ia sample is dominate by 
true protostars. We could then be witnessing the 
onset of X-ray activity in very young YSOs. 

Class O-I objects are affected by larger absorp- 
tions with respect to the other classes as we will 
confirm in Sect. 15. 3] using the X-ray derived Nh 
values. Can this fact results in a systematic un- 
derestimation of Lx values for Class O-I stars? 
Although we do apply absorption corrections to 
the observed X-ray fluxes when deriving luminosi- 
ties both in the 0.5-8.0 keV and 2.0-8.0 keV bands, 
these are corrections to the observed X-ray spectra 
and emission from relatively cool plasma would re- 
main unobserved, and therefore unaccounted for, 
in heavily observed sources. 

In order to verify that the X-ray luminosities of 
heavily absorbed Class O-I objects are accurately 
corrected for absorption, we assumed that these 
sources have X-ray spectra similar to Class II and 
III objects, and computed through extensive sim- 
ulations the effect of varying the absorbing hydro- 
gen column density on the determination of X-ray 
luminosities. 

For this purpose we considered the best fit spec- 
tra of a subsample of ^200 Class II and III ob- 
jects with small absorption (LogNn < 21.5cm~^), 
consisting of either one or two thermal compo- 
nents. For each of these spectra, using XPSEC 
and the fakeit command, we computed 10 simu- 
lated spectra for each of five different Nh values in 
the range lO^^cm^^ <Nh < lO^^cm"^. We then 
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analyzed the simulated sp ectra with a procedure 
analogous to that used by iGetman et all (|2005bh 
for the analysis of the observed spectra. Each 
spectrum was rebinned so to have a minimum of 
15 photon per channel and was then fit with both 
an isothermal and a two component model. After 
adopting the simplest of the two models that gave 
a statistically acceptable fit, we computed the un- 
absorbed flux in the 0.5-8.0 keV and 2.0-8.0 keV 
bands. We then computed the absorption cor- 
rected XLFs for each of the five absorption values. 
Most of the simulated sources did not have enough 
photons when absorbed by Nh = lO^'^cm"^. In 
order to avoid selection effects, since we are not 
accounting for upper limits, we could only com- 
pare the XLFs for the 115 sources whose spectra 
were bright enough to be analyzed even when con- 
sidering the highest absorption value. The result 
is that the 5 XLFs for the 0.5-8.0 keV luminosities 
have median within ~Q.ldex but tend to become 
wider as the Nh is increased. The five XLFs for 
2.0-8.0 keV luminosities are instead very similar 
in shape and with medians within 0.03 dex, ex- 
cept for the one relative to the highest absorption 
values which is brighter hy ~0.1 dex. We conclude 
that, in the assumption that Class O-I sources have 
the same spectra of Class ILIII ones, the X-ray lu- 
minosities we derive are not underestimated. 

We cannot exclude that some of the Class O-Ib 
stars are actually heavily absorbed Class II stars 
(see Fig. [3]). However, even in this case, our result 
would plausibly indicate a significant evolution of 
X-ray activity in the early stages of YSO forma- 
tion. 

Due to various IR sensitivity limitations at very 
high obscurations, the Class O-I sampling might be 
restricted to higher masses (with stronger IR bolo- 
metric luminosities), and through the Lx-mass 
correlation, would have higher Lx values introduc- 
ing a bias. But we find that these heavily obscured 
stars have lower, not higher, Lx values, further 
supporting our result. 

Another bias may be introduced by the require- 
ment that all stars considered have IRAC detec- 
tions in one or more bands. Therefore, it is pos- 
sible that lower-mass (i.e. bolometrically weaker) 
Class III stars have been systematically omitted 
from the sample, leaving behind too many higher 
mass, i.e. stronger Lx stars. But this is not rele- 
vant as the XLF comparisons in Figs.[8]and[9]are 



binned by mass. We might have few 0.1-0.5 
Class III stars, but those shown should be unbi- 
ased in Lx. 

Another possible worry is that flares could dom- 
inate the flux from some of our objects, and our 
X-ray spectra could pertain to the flaring state 
rather than the quiescent state. We have verified 
that this is not the case by u sing the technique 
described in IWolk et al.l (|2005( ) to recompute X- 
ray spectra and fluxes for our stars using photons 
detected only during non-flaring intervals. Fits to 
these new spectra yield XLFs that do not differ 
significantly from those shown in Fig.[8l We have 
also verified that using these new spectra the re- 
sults discussed in the next section do not change 
significantly. 

We note that our results are consistent with 
the non-detection of class sources of Serpens 
( Giardino et al.ll2007 ). In fact, their limiting sensi- 
tivity, (4.0 X 10^^ erg/s in the hard band) estimated 
assuming temperature and absorption similar to 
those we found for our X-ray detected Class 0- 
I sources, corresponds to the highest tenth per- 
centile of the XLF of our Class 0-Ia stars. Con- 
sidering the small size of the Serpens Class star 
sample, their non-detection is not surprising. 

5.3. Gas and dust absorption 

Figure[n] shows the Spitzer/IRAC [3.6]-[4.5] 
color vs. Chandra X-ray median energy of all se- 
lected sources. The typical uncertainty in the me- 
dian energy has been estimated us ing the results 
of the MARX simulations made in iGetman et al.l 
( 20061 ). where the median energy errors were com- 
puted as a function of the source net counts. We 
used the value of 0.3 keV, corresponding to objects 
with MedE>2keV and with about 30 net counts 
which is the minimum of net counts we have for 
Class O-I objects. The typical uncertainty in the 
IRAC color has been estimated assuming a 5% 
uncertainty for the IRAC magnitudes. 

2007h . high X- 



As discussed in IGetman et al.l 
ray median energy (MedE> 1.7 keV) can be con- 
sidered as an indicator of absorbing column den- 
sity that arises from gas, primarily metallic ele- 
ments (e.g. N, O, Ne, Mg, Si, Ar, Fe, etc.). On 
the other hand, if the object is known to be a 
protostar, the [3.6]-[4.5] color is primarily a mea- 
sure of the density of grains in the local infalling 
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Fig. 11.— Spitzer/mkC [3.6]-[4.5] color vs. 
Chandra X-ray median energy of all selected 
sources. The different classes ONC objects are in- 
dicated using the same symbols of Fig.[5l while 
filled squares are t he values of the Y SO in 
IC 1396N taken from lGetman et all (|2007[) . Typ- 
ical error bars in the median energy and in the 



IRAC colors are also shown. 

envelope. Thus, the cor relation between thes e two 
quantities, first found in lGetman et al.l (j2007l ). has 
been interpreted by these authors as evidence of 
absorption from the protostellar envelope rather 
than from the ambient molecular cloud material. 
Although our data do not cover the whole X-ray 
median energy, they show a trend consistent with 
that found using the IC 1396N data, also shown 
in Fig. 111! Nevertheless, the interpretation of our 
data is not as clear as it was in IC1396N, since 
the ONC background is quite complex while the 
maximum line-of-sight column is Ay ^ 10 mag. 

Fig. [12] compares the absorbing column density 
Nh for the different populations. In the left panel, 
logNn is plotted for the X-ray detected objects 
as a function of the unabsorbed X-ray luminosity 
Lt.c- The right panel shows, for each class, the me- 
dian values of the absorbing column density distri- 



butions, values between the first and third quar- 
tiles and values between the minimum and the 
maximum of the distributions of the logNn val- 
ues. Note that only objects with errors in logNn 
smaller than 0.5 have been considered. How- 
ever, the median values of the distribution do not 
change significantly if all data are considered. 

A clear trend of increasing absorbing column 
density from Class III to Class O-I can be seen both 
in the left and right panels. This trend is expected 
since Class O-I objects are characterized by a large 
amount of infalling material and are also typi- 
cally embedded in the molecular cloud; a smaller 
amount of obscuring material is instead present for 
Class II objects that have shed their envelope and 
formed a circumstellar disk, while the extinction 
suffered from Class III objects is due to interstellar 
material only (both within the cloud and on the 
line of sight to the cloud). 

We note that the estimate of the Nh values of 
Class O-I stars is based primarily on the hard com- 
ponent of the X-ray spectrum. That is, at the 
temperatures of the protostars (kTav ^ 2keV) 
and with log Nh ~ 23cm~^, the observed flux 
is almost entirely due to the hard component of 
the X-ray spectrum. At the same temperatures, if 
the Nh was smaller, the total flux would be dis- 
tributed between the soft and hard component and 
the shape of the spectrum would be different. Our 
technique should thus be the most accurate for ob- 
jects with high Nh values - thus yielding particu- 
larly good Nh values for the, on average, heavily 
embedded protostars. Therefore, we argue that 
the low XLFs for protostars we have derived (see 
Figure 9) are unlikely to be due to errors in our 
estimates of the X-ray absorption towards these 
objects, as confirmed from the analysis of the hard 
XLFs and the simulated spectra of Class II and III 
(see Sect.lO]). 

5.4. Evolution of YSO X-ray spectra and 
temporal characteristics 

Fig. [13] shows the results of the comparison of 
the average plasma temperature fcTav amongst the 
different classes. For X-ray spectra fit with two 
isothermal components, fcTav was computed as the 
mean of the soft and hard fcT spectral components, 
weighted by their respective emission measures; 
for spectra fit with a single component fcTav=fcT. 
The fcTav values for the X-ray detected objects 
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Fig. 12. — Left panel shows the logNn values as a function of the unabsorbed X-ray luminosity Lt_c for the 
X-ray detected objects. Symbols are as in Fig.O Right panel shows, for each class, the median values of 
the absorbing column density distributions indicated by the big black dots. Solid bars indicate the values 
of logNn included between the first and third quartiles, while the dashed lines (partially hidden by solid 
bars) indicate the values included between the minimum and the maximum of the distributions of the log Njj 
values. Note that only objects with errors in logNjj smaller than 0.5 were considered. 



22 



of the diflFerent populations are plotted as a func- 
tion of the unabsorbed X-ray luminosity T^t,c- Note 
that only objects with errors on /cTav smaller than 
1 keV have been considered. However the distribu- 
tions do not change significantly if all data are con- 
sidered. The bottom-right panel of Fig. [T3l shows, 
for each class, the median values of the average 
plasma temperatures, values between the first and 
third quartiles and values between the minimum 
and the maximum of the distributions of the fcTav 
values. No trend in plasma energy is found as a 
function of evolutionary stage from Class 0-Ia to 
Class III systems. 

A trend of increasing of the average plasma 
temperature with increasing unabsorbed X-ray lu- 
minosity is clearly evident for all the samples, as 
has often been r eported for other ty pes of coro- 
nal sources (e.g. ISchmitt et al. 19901 ). A similar 
result is obtained by considering only the sub- 
sample of stars with X-ray spectra fit with two 
temperature components: high luminosity stars 
tend to have high values of the ratio between 
the high and low temperature emission measures 
(EM2/EM1 ^ 2.5), suggesting that the hot plasma 
component dominates over the cold one for the 
high luminosity stars. No significant difference be- 
tween the classes is observed in Fig.fTSl 

Finally we checked if ONC stars of different 
classes show different X-ray temporal behavior. 
To do this, we used the results of the nonparamet- 
ric one-samp le Kolmogorov - Smirno v (K-S) test 
performed bv ICetman et aL I (|2005d) . The results 
of the K-S test, Pks, strongly depend on the 
source count-statistics and therefore on their count 
rate. For this reason, we compared the probabil- 
ities of variability Pvar=(l-PKs) for the different 
classes of stars with count rates in definite ranges. 
The results are given in Tabled where, for each 
class and for a given range of Count Rate (CR), we 
report the number of stars with Pvar > 99% with 
respect to the total number of stars of each class in 
the corresponding range of count rates. The per- 
centage of variable stars is also indicated within 
brackets. Even if the statistics for each sample is 
quite poor, the percentages indicate that the dif- 
ferent classes have very similar fractions of variable 
stars. In practice, when observed with high statis- 
tics, all the sources are variable, independently on 
their class. 



6. CLASS O-I: COUNTERPARTS AT OTHER 
WAVELENGTHS 

We list in TablcfTil the counterparts of the 
Class O-I objects identified in this work that have 
been detected in X-rays ( the COUP survey). 



JHKL bands l|Lada et all l2000t iMuench et al] 



2OO2I hereafter LMHOO and MLL02, re spectively), 
10 and 20 /zm ([Robberto et al.J 2005a, hereafter 
RBP5) and in the 11.7/im (jSmitheLalJ l2005l 



hereafter SBS05). None of our candidate Class O-I 
sources are in the li st of 1.3 cm radio sources of 
Zapata et af] (|2004bl ). 



Th e JHKL photometric survey from lLada et al 



(|2000l ) included 391 stars in the central 6!55x6;55 
of the Trapezium cluster. By selecting stars with 
K — L >1.5, they obtained a list of 78 candidate 
protostars. We correlated this list with our cata- 
log using a ma tching radius of 1''2. As expected, 
the Lada et al. I (I2OOQ) sample, selected on the ba- 
sis of the K — L color, is larger than our sample 
for which the selection requires at least two colors 
computed using the K and the four IRAC magni- 
tu des. For thi s reason, we find that only 4 of the 
78 Lada et al. ( 2000f ) candidate protostars are in 
common with our Class 0-Ia stars and 9 with our 
Class 0- lb stars. Among the remaining 65 objects, 
16 cannot be Class O-I objects since they have an 
optical counterpart in our catalog and 27 have 
[3. 6] - [4. 5] < 0.7 and are therefore consistent with 
star-disk systems (Cl ass II) rathe r than protostars; 
for the remaining 22 iLada et all (|2000l ) candidate 
protostars, the available magnitudes do not satisfy 
our selection criteria. 

Among our Class O-I objects, only 5 (Class 0- 
Ib) have been detected at 10 and 20 ^m by 
Robberto et al. I (|2005al) : of these 5, 2 are in the 
lis t of the embedded objects previously studied 
(|2004L 



Smith et al 



hereafter SBS04) in the 
OMC-IS region at 8.8, 11.7 and 18.8 /xm. 

We i nterpolated the flux conversion coefficients 



given in lLada et al.l (|2006i) to 10 and 20 /im and we 



computed the corresponding fluxes. SEDs with all 
the available fluxes for these 5 objects are shown 
in Fig.IH 

In order to understand the properties of our 
sample of Class O-I objects, we also c onsidered the 
compac t rad io sources detected b y IZapata et al 



(|2004bl) and IZapata et al. I dlooii) usmg, respec- 
tively, 1.3cm (in a region of 30"x30"around 



23 




Fig. 13. — kTav values as a function of the unabsorbed X-ray luminosity Jjt.c for the X-ray detected objects, 
plotted separately for Class O-I, Class II and Class III stars, in the first three panels. Symbols are as in 
Fig. [51 The bottom-right panel shows, for each class, the median values of the average plasma temperature 
distributions indicated by the big black dots. Solid bars indicate values of the temperature distributions 
between the first and third quartiles, while dashed lines (partially hidden by solid bars) indicate values 
between the minimum and the maximum of the distributions of the kTav values. Note that only objects 
with errors in kTav smaller than 1 keV were considered. 
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Table 4: X-ray variable sources of different evolutionary class stars in the ONC with Pks > 99%. 



log (CR) 


Class 0/Ia 


Class 0/Ib 


Class II 


Class III 


[-5,-4] 


2/ 4 ( 50%) 


0/ 2 ( 0%) 


1/ 4 ( 25%) 


0/ 2 ( 0%) 


[-4,-3 ] 


5/ 6 ( 83%) 


6/ 9 ( 66%) 


15/28 ( 53%) 


19/30 ( 63%) 


[-3,-2 ] 




7/ 7 (100%) 


38/38 (100%) 


46/50 ( 92%) 



[3.6] [4.5] [5. a] [8.0] 




[3.6] [4.5] [5.8] [a.O] 
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Fig. 14. — SED of Class O-I candidate objects in the lRobberto et al.l (|2005al ) list with photometry at 10 and/or 
20 fim. Dots indicate fluxes computed used the magnitudes fro m our catalog, plus symbols correspond to the 
magnitudes at 10 and 20 fim from the Robberto et al.l (j2005al ) catalog and X correspond to the magnitudes 
at 8.8, 11.7 and 18.8 /im from the Smith et al. ('2004') catalog. Multiple values at 10 /im are relative to the 
two different nights of observation in .Robberto et al.. (,2005a) . 
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OMC-IS) and 3.6 cm observations (in a region 
of 4' X 4' around the Trapezium cluster) performed 
with the Very Large Array. With a matching ra- 
dius of 1", we retrieved in our catalog 7 of the 11 
radio sources at 1.3 cm and 67 of the 77 compact 
3.6 cm sources. 

None of the IZapata et all (|2004al lbh matched 
sources is in our list of Class O-I objects. The 
absence of 1 cm radio counterparts in our proto- 
star sample might imply that the radio sources are 
mainly magnetically active non-thermal Class III 
stars rather than thermal protostellar outflows. 

Note that the north-east region of the COUP 
FOV, where most of our candidate Class O-I are 
concentrated, is not included in the FOV of the 
literature sources mentioned in this section. For 
this reason, the SEDs of Class O-I objects in this re- 
gion cannot yet be extended to wavelengths longer 
than 8 /im, which is crucial to confirm (or refute) 
their protostellar nature. 

Finally, our list of candidate Class O-I proto- 
stars includes only two objects (COU P sources 420 
and 48 4) of the list of stars studied in lGrosso et al. 
(|2005f ). located at the bottom edge of the OMC- 
IS region, i.e., far from the density peak and its 
associated nebular emission. We know that the 
absence of IRAC protostars near OMC-ls and in 
the BN/KL region is due to nebular contamination 
and crowding, not to the real absence of protostars 
in these two high density regions. In addition, we 
cannot classify the COUP sources 5 82, 594, 615 
633, 6 41, 659 and 667, indicated in iGrosso et al 
as embedded YSOs in a subcluster, since 



they do not have any counterpart at the IRAC 
bands. 

7. SUMMARY AND CONCLUSIONS 

We have selected a sample of 45 candidate 
Class O-I objects in the ONC field observed in the 
X-rays by COUP using new deep JHK, [3.6], 
[4.5], [5.8] and [8.0] photometry obtained at the 
ISPI@4m CTIO telescope and from Spitzer IRAC 
observations. New deep optical observations taken 
with the WFPC2 camera of HST and the WFI 
camera of ESO 2.2 m were also used to reject can- 
didate Class O-I objects with optical counterparts. 

We distinguish between candidate protostars 
with SEDs increasing monotonically from K up to 
8 fim, indicated as Class 0-Ia, and candidate pro- 



tostars with SED increasing from K up to 4.5 /xm 
and decreasing at longer wavelengths, indicated 
as Class 0-Ib, since the different slope of the SED 
at wavelengths larger than 4.5 /xm could be an in- 
dicator of a different density in the circumstellar 
material. 

Out of the 23 Class 0-Ia stars, 10 have been de- 
tected in the X-rays with COUP, while out the 
22 Class 0-Ib stars, 18 have been detected; we are 
thus in the position to study, for the first time, the 
X-ray properties of a significant and homogeneous 
sample of candidate protostars. 

Samples of bona fide more evolved Class II and 
III stars in the same region have also been selected 
in order to study the time and mass evolution of 
X-ray activity. 

Our principal result is that Class 0-Ia stars 
are significantly less luminous in X-rays, both 
in the total and hard bands, than more evolved 
Class II stars with mass larger than 0.5 Mq, while 
Class 0-Ib stars show X-ray luminosities similar to 
Class II stars; finally, our data confirm previous 



determinations (|Neuhauseiill997t iFlaccomio et al . 
2003bl : IPreibisch et all l2005h that Class III stars 



with masses >0.3Mq are more X-ray luminous 
than Class II objects, even if the result is quite 
marginal with our data. 

Our result that Class 0-Ia objects are less lu- 
minous in X-rays than objects of other classes, 
support the hypothesis that the onset of X-ray 
emission occurs at a very early stage, when the 
objects show an increasing SED at least up to 
8.0 iim. The lack of a detailed theory of proto- 
stellar evolution does not allow us to interpret 
the difference in the evolutionary stage between 
Class 0-Ia and 0-Ib protostars. Although Class 0- 
Ib objects are selected as candidate protostars, our 
data suggest that this sample could be dominated 
by reddened Class II stars. If we consider the X - 
ray limiting sensitivity in iGiardino et al.l (j2007l ). 



we deduce that our results are in agreement with 
the non-detection of Class protostars in Serpens. 

The X-ray spectral properties of Class 0-Ia and 
0-Ib are similar to those of the more evolved 
Class II and III objects except for a larger absorp- 
tion likely due to the enhanced local density sur- 
rounding protostellar objects. The different evo- 
lutionary classes also show similar X-ray variabil- 
ity characteristics and the protostellar variabil- 
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ity is indistinguishable from T Tauri variability. 
This result supports the XEST determination that 
variability of Class II and III stars are the same 
(fStelzer et al.ll2007h and the COUP determination 
that T Tauri M star s have the same variab ility as 
more massive stars (jCaramazza et al.ll2007r ). Also 
related is the COUP result that YSOs with high 
X-ray luminosities have the same first ionization 
potential-rel a.ted elemental abu ndances as older 
flaring stars (jMaggio et al.ll2007l ). 

All of these results support the general theme 
that there is only one mechanism in YSO X- 
ray production at energies >lkeV, and that this 
mechanism arises from stellar magnetic activity 
which does not depend on ac cretion or disk s . The 
various receiit arguments dKastner et al. 2002t 



Schmitt e t ar2005:'Gi idel et al ."2007':'Argiroffi et al 
20071 : iHuc nemocrd er et al.l 2007, : .Grosso et al 
20071 ) that accretion accounts for some X-ray emis- 



sion and line ratios, probably only refers to soft 
components which are not accessible to our study 
(|Feigelson et al.ll2007l) . 
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Table 5: NIR photometry of Class 0/Ia candidates. 
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Table 6: NIR photometry of Class 0/Ib candidates. 
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11.515± 0.007 10.718± 0.004 10.300± 0.030 










25 




83. 


.7545867 


-5, 


.3937160 












14, 


.382± 


0, 


.011 


11.811± 0.005 11.018± 0.005 10.743± 0.017 










26 


167 


83, 


,7570887 


-5, 


.4201451 












15, 


.427± 


0, 


.020 


12.991± 0.014 12.173± 0.013 11.450± 0.036 










27 


209 


83, 


.7672732 


-5, 


.3755507 




15. 


.475± 


0, 


.013 


13, 


.013± 


0, 


.016 


10.755± 0.004 


9.768± 0.003 


9.082± 


0.007 


8, 


,151± 


0, 


.018 


28 


407 


83, 


.7957359 


-5, 


.2559454 16.583± 0.043 12. 


.803± 


0, 


.037 


10, 


.836± 


0, 


.032 


9.232± 0.001 


8.464± 0.004 


8.104± 


0.004 


7. 


.707± 


0, 


.033 


29 


420 


83, 


.7970836 


-5, 


.4106294 




16. 


.107± 


0, 


.023 


12, 


.677± 


0, 


.013 


9.662± 0.008 


8.171± 0.004 


7.447± 


0.014 










30 


448 


83, 


.7993539 


-5, 


.3500809 




12, 


.617± 


0, 


.024 


10, 


.011± 


0, 


.011 


7.644± 0.002 


6.807± 0.001 


6.127± 


0.005 


5, 


.326± 


0, 


.012 


31 


472 


83, 


.8007853 


-5, 


.4156826 




12, 


.477± 


0, 


.009 


10, 


.698± 


0, 


.012 


8.904± 0.003 


8.157± 0.004 


7.487± 


0.013 


6, 


.636± 


0, 


.027 


32 




83, 


.8058429 


-5, 


.3914514 




















8.621± 0.009 


7.010± 0.007 


5.956± 


0.008 


5, 


.096± 


0, 


.035 


33 


548 


83, 


.8064706 


-5, 


.2960711 




14, 


.648± 


0, 


.009 


11, 


.908± 


0, 


.012 


9.992± 0.002 


9.181± 0.002 


8.590± 


0.010 


7, 


.577± 


0, 


.019 


34 


570 


83, 


.8071662 


-5, 


.3599724 




12, 


.897± 


0, 


.023 


10, 


.624± 


0, 


.012 


8.404± 0.011 


7.445± 0.010 


6.960± 


0.032 










35 




83, 


.8074223 


-5, 


.3944417 




















7.047± 0.005 


5.170± 0.006 


3.888± 


0.004 


2, 


.994± 


0, 


.038 


36 


595 


83, 


.8079975 


-5, 


.4501972 




15, 


.196± 


0, 


.014 


11, 


.901± 


0, 


.021 


9.232± 0.004 


8.236± 0.001 


7.752± 


0.012 


7, 


.369± 


0, 


.032 


37 




83, 


.8195534 


-5, 


.3125430 












13, 


.387± 


0, 


.046 


9.792± 0.008 


8.530± 0.002 


7.756± 


0.018 


7, 


.093± 


0, 


.042 


38 


1030 


83, 


.8307100 


-5, 


.4283582 17.377± 0.083 13. 


.356± 


0, 


.022 


10, 


.409± 


0, 


.016 


7.654± 0.002 


6.779± 0.001 


6.034± 


0.009 


5, 


.266± 


0, 


.021 


39 


1054 


83, 


.8325174 


-5, 


.2597844 




16, 


.173± 


0, 


.059 


12, 


.783± 


0, 


.030 


9.339± 0.002 


8.100± 0.004 


7.200± 


0.002 


6, 


.411± 


0, 


.016 


40 


1094 


83, 


.8346095 


-5, 


.3326871 




13, 


.958± 


0, 


.014 


11, 


.398± 


0, 


.008 


8.927± 0.005 


8.073± 0.004 


7.824± 


0.032 










41 


1170 


83, 


.8408440 


-5, 


.3354363 




14, 


.985± 


0, 


.016 


12, 


.799± 


0, 


.017 


10.155± 0.006 


9.281± 0.005 


9.164± 


0.110 










42 


1188 


83, 


.8420775 


-5, 


.3160413 




16, 


.641± 


0, 


.021 


13, 


.638± 


0, 


.014 


10.825± 0.016 


9.652± 0.006 


8.930± 


0.047 










43 


1364 


83, 


.8604067 


-5, 


.3218452 




15, 


.505± 


0, 


.018 


12, 


.678± 


0, 


.010 


10.317± 0.008 


9.364± 0.004 


8.923± 


0.032 










44 


1480 


83, 


.8822676 


-5, 


.3556764 




14, 


.409± 


0, 


.018 


11, 


.755± 


0, 


.016 


9.401± 0.004 


8.554± 0.002 


8.868± 


0.032 










45 


1565 


83, 


.9236385 


-5, 


.3374070 












15, 


.703± 


0, 


.019 


12.773± 0.013 11.751± 0.008 


11.234± 0.035 











Table 7: NIR photometry of Class II candidates. 


Seq. COUP RA(2000) Dec(2000) J 
Num. ID [deg] [deg] 


H 


K 


[3.6] 


[4.5] 


[5.8] 


[8.0] 



46 
47 
48 
49 
50 
51 
52 
53 
54 
55 
56 
57 
58 
59 
60 



12 
13 
21 

29 



37 
40 
44 

54 
55 
58 
62 
66 



83.6700479 
83.6706135 
83.6853649 
83.6939237 
83.6948357 
83.6949818 
83.6996888 
83.7003928 
83.7044071 
83.7082042 
83.7101265 
83.7104375 
83.7113833 
83.7148246 
83.7167237 



-5.4440551 
-5.4332919 
-5.4106835 
-5.3903819 
-5.4346142 
-5.3580154 
-5.3944798 
-5.3772892 
-5.4406884 
-5.3123055 
-5.3389923 
-5.3930997 
-5.4502169 
-5.4201915 
-5.4119330 



12.550± 0.024 11.652± 0.043 
13.445± 0.018 12.312± 0.022 
13.558± 0.013 12.778± 0.010 
12.720± 0.013 11.771± 0.013 
14.469± 0.024 12.810± 0.025 
13.446± 0.013 12.466± 0.013 



13.878± 0.013 
13.248± 0.026 
8.795± 0.041 
12.107± 0.030 
12.492± 0.024 
11.527± 0.030 



13.138± 0.012 
12.420± 0.022 
8.383± 0.048 
11.183± 0.015 
11.626± 0.029 
10.750± 0.034 



12.304± 0.037 11.256± 0.036 



11.089± 0.029 
11.533± 0.019 
12.406± 0.013 
11.112± 0.013 
11.842± 0.023 
11.948± 0.012 

12.622± 0.010 
12.035± 0.017 
7.687± 0.085 
10.717± 0.016 
11.202± 0.019 
10.445± 0.022 
9.656± 0.020 
10.730± 0.035 



10.416± 0.005 
10.571± 0.005 
11.811± 0.007 
10.240± 0.002 
10.443± 0.002 
11.151± 0.009 
12.039± 0.006 
12.052± 0.006 
11.302± 0.004 
6.022± 0.001 
9.871± 0.002 
10.872± 0.003 
9.557± 0.001 
8.562± 0.001 
9.748± 0.002 



10.118± 0.006 
10.201± 0.006 
11.353± 0.010 
9.668± 0.004 
9.940± 0.002 
10.750± 0.007 
11.690± 0.007 
11.749± 0.010 
10.854± 0.004 
5.520± 0.001 
9.532± 0.003 
10.517± 0.004 
9.214± 0.002 
8.092± 0.001 
9.289± 0.001 



9.754± 0.014 
9.904± 0.014 
11.119± 0.032 
9.113± 0.009 
9.608± 0.007 
10.291± 0.037 
11.278± 0.027 
11.367± 0.027 
10.541± 0.014 
5.129± 0.002 
9.279± 0.010 
10.190± 0.015 
8.926± 0.004 
7.725± 0.002 
8.944± 0.005 



9.001± 0.016 
9.172± 0.017 
10.430± 0.107 
8.115± 0.027 
8.951± 0.011 
9.525± 0.060 
10.356± 0.046 
10.517± 0.079 
9.691± 0.022 
4.470± 0.001 
8.605± 0.023 
9.353± 0.029 
8.151± 0.006 
6.816± 0.002 
8.198± 0.006 



Table 8: NIR photometry of Class III candidates. 


Seq. COUP RA(2000) Dec(2000) J 
Num. ID [deg] [deg] 


H 


K 


[3.6] 


[4.5] 


[5.8] 


[8.0] 



194 1 83.6227896 

195 83.6230879 

196 83.6400806 

197 6 83.6592044 

198 83.6594147 

199 83.6614465 

200 83.6650038 

201 7 83.6656279 

202 9 83.6661159 

203 10 83.6666439 

204 83.6692738 

205 14 83.6734241 

206 16 83.6743312 

207 83.6760625 

208 83.6761015 



5.3936165 14.027± 0.011 
5.3954869 12.020± 0.029 
•5.3860815 11.285± 0.027 
■5.4065827 12.835± 0.012 
•5.4007576 15.957± 0.014 
•5.3569493 14.159± 0.016 
5.4217564 18.447± 0.044 
5.4071390 8.751± 0.034 
•5.4449746 10.138± 0.025 
5.4343872 13.298± 0.018 
5.3876238 16.189± 0.013 
•5.3992351 13.280± 0.012 
5.3637129 14.047± 0.013 
5.3758147 16.756± 0.017 
•5.4214457 15.926± 0.018 



13.316± 0.013 
10.659± 0.017 
10.754± 0.015 
11.990± 0.013 
14.627± 0.014 
13.023± 0.016 
16.860± 0.039 
8.408± 0.071 
9.908± 0.017 
12.589± 0.017 
15.246± 0.014 
12.558± 0.013 
13.394± 0.015 
15.252± 0.017 
14.885± 0.020 



13.092± 0.011 
10.134± 0.015 
10.648± 0.014 
11.698± 0.017 
14.185± 0.017 
12.588± 0.021 
16.321± 0.023 
8.130± 0.079 
9.355± 0.016 
12.240± 0.015 
14.938± 0.027 
12.291± 0.013 
13.009± 0.013 
14.567± 0.019 
14.475± 0.019 



12.757± 0.019 
9.780± 0.004 
10.743± 0.004 
11.544± 0.007 
14.127± 0.035 
12.409± 0.014 
15.897± 0.100 
7.890± 0.002 
9.417± 0.003 
11.951± 0.009 
14.723± 0.051 
11.952± 0.012 
12.587± 0.025 
14.280± 0.074 
14.164± 0.029 



12.705± 
9.642± 
10.776± 
11.417± 
14.129± 
12.446± 
16.088± 
7.895± 
9.413± 
11.869± 
14.558± 
11.852± 
12.510± 
14.171± 
14.068± 



0.024 
0.004 
0.007 
0.010 
0.062 
0.020 
0.182 
0.003 
0.004 
0.013 
0.069 
0.016 
0.023 
0.101 
0.041 



9.625± 0.014 9.586± 0.026 
10.794± 0.020 10.685± 0.069 
11.525± 0.033 



12.528± 0.141 

7.838± 0.005 
9.400± 0.010 
11.863± 0.046 

11.862± 0.124 



7.772± 0.007 
9.375± 0.019 



Table 9: X-ray photometry of the Class O/I candidates taken from iGetman et al.l (|2005br ). A description of 
the columns is given in Section [5T2l 



Seq. COUP NetCts PSF Exp. log Nh kTl log EMI log Lh^c log Lt,c 



Num. 


iU 




Frac 


(ks) 




cm 


J 




[Key ) 




cm 




(erg 


s ) 


(ergs 1) 


o 



A 

Zli 


i (U 





86 


762 


20 


22 


56± 


U 


Oo 


1 c nn_L 1 t: nn 
i5.UU± io.UU 


5z 


oo 1 

2o± 





06 


on 
z9 


oo 
Z» 


29 


42 


1 


484 


54 





86 


783 


40 


22 


54± 


U 


1 A 

14 




51 


9d± 





31 


oo 
ZS 


<0l 


29 


00 


9 


d9d 


/ 1 





87 


825 


80 


22 


i6± 


U 


oo 
Z<5 


y.L)D± io.UL) 


51 


i9± 





31 


oo 
ZS 


1 n 
19 


28 


38 


ID 


luZ 


(^9 





86 


801 


10 


oo 
22 


31± 





1 1 


J.z5± i. ( 1 


51 


o O 1 

<3D± 





27 


oo 
ZS 


63 


28 


93 


ii 


860 


IOC 
1<35 





87 


792 


20 


oo 
22 


60± 





1 o 


5. (o± 4.uy 


5z 


o4± 





20 


on 
z9 


2( 


29 


49 


1 o 


o c n 
(359 


4oi 





87 


751 


60 


oo 
22 


1 

1 1± 





05 




53 


oo 1 

zoi 





12 


on 
z9 


10 


30 


24 


16 


1048 


162 





87 


795 


80 


22 


68± 





08 


2.98± 0.84 


52 


51± 





19 


29 


24 


29 


56 


17 


1115 


27 





86 


794 


00 


22 


83± 





19 


1.13± 0.70 


52 


43± 


1 


00 


28 


56 


29 


37 


19 


1197 


15 





87 


799 


30 




























21 


1321 


184 





87 


714 


40 


23 


11± 





09 


4.51± 2.98 


52 


82± 





29 


29 


69 


29 


94 


24 


103 


58 





87 


783 


40 


22 


63± 





13 


0.95± 0.35 


52 


64± 


1 


00 


28 


61 


29 


60 


26 


167 


88 





86 


739 


20 


22 


73± 





10 


1.16± 0.24 


52 


79± 


1 


00 


28 


94 


29 


72 


27 


209 


97 





87 


824 


10 


22 


37± 





13 


4.86± 3.60 


51 


90± 





22 


28 


80 


29 


03 


28 


407 


479 





87 


737 


40 


22 


63± 





07 


3.61± 0.99 


52 


91± 





16 


29 


71 


30 


00 


29 


420 


65 





86 


838 


20 


22 


29± 





13 


3.45± 2.65 


51 


70± 





26 


28 


48 


28 


78 


30 


448 


2249 





86 


820 


50 


22 


52± 





02 


2.24± 0.12 


53 


64± 





05 


30 


24 


30 


64 


31 


472 


501 





86 


456 


20 


22 


55± 





05 


2.80± 0.38 


53 


18± 





10 


29 


88 


30 


21 


33 


548 


2315 





87 


774 


50 


22 


54± 





02 


3.06± 0.25 


53 


59± 





04 


30 


33 


30 


65 


34 


570 


4070 





87 


829 


40 


22 


64± 





02 


3.35± 0.22 


53 


86± 





03 


30 


63 


30 


93 


36 


595 


624 





87 


779 


80 


22 


69± 





04 


2.31± 0.18 


53 


23± 





10 


29 


84 


30 


23 


38 


1030 


3616 





87 


749 


80 


22 


98± 





02 


15.00± 1.00 


53 


89± 





03 


30 


88 


31 


01 


39 


1054 


1630 





88 


707 


30 


22 


64± 





02 


2.59± 0.25 


53 


62± 





07 


30 


28 


30 


64 


40 


1094 


214 





87 


801 


10 


22 


61± 





07 


2.79± 0.55 


52 


60± 





17 


29 


30 


29 


64 


41 


1170 


150 





87 


799 


30 


22 


61± 





08 


1.03± 0.22 


52 


96± 


1 


00 


29 


02 


29 


91 


42 


1188 


232 





87 


686 


10 


22 


73± 





04 


1.70± 0.27 


53 


10± 





16 


29 


54 


30 


05 


43 


1364 


501 





88 


779 


80 


22 


70± 





03 


4.30± 1.02 


52 


94± 





11 


29 


80 


30 


05 


44 


1480 


914 





86 


684 


40 


22 


69± 





03 


2.10± 0.18 


53 


51± 





08 


30 


07 


30 


50 


45 


1565 


385 





88 


746 


20 


22 


52± 





05 


2.50± 0.40 


52 


85± 





12 


29 


50 


29 


86 



34 



Table 10: Upper Limits of the X-ray undetected Class O/I candidates. A description of the columns is given 
in Sectionl5.2l 



Seq. 


RA (2000) 


Dec (2000) 


Cts 


Cts 


Exp. 


log 


T 1 


log 


^h.c 


log 




log 




log 


Lie 


log 


Lie 


Num. 




(deg) 




(deg) 


hard 


tot. 


(ks) 


(erg 


s-1) 


(erg 


s 1) (erg 


s-1) 


(erg 


s 1) (erg 


s-i) 


(erg 


s-i) 


1 


83, 


,70931200 


-5, 


.36415800 


26.38 


32.58 


710.29 


28, 


.44 


28, 


.27 


28 


.76 


28, 


.84 


29, 


.44 


29 


.65 


2 


83, 


,73136100 


-5, 


.29824800 


69.32 


86.29 


739.17 


28, 


.84 


28, 


.67 


29 


.16 


29, 


.25 


29, 


.85 


30 


.05 


4 


83, 


,79406000 


-5, 


.58508600 


661.04 


3605.31 


681.67 


29, 


.85 


29, 


.69 


30 


.18 


30 


.90 


31 


.50 


31 


.71 


5 


83, 


,79568709 


-5, 


.36277907 


18.18 


31.25 


817.28 


28, 


.21 


28, 


.05 


28 


.54 


28, 


.76 


29, 


.36 


29, 


.57 


6 


83, 


,80009500 


-5, 


.31133400 


21.00 


26.16 


781.40 


28, 


.29 


28, 


.13 


28 


.62 


28, 


.70 


29, 


.30 


29, 


.51 


8 


83, 


,81112700 


-5, 


.28122600 


30.07 


64.72 


632.86 


28, 


.54 


28, 


.38 


28, 


.87 


29, 


.19 


29 


.79 


30 


.00 


13 


83, 


,82558451 


-5, 


.34145849 


16.28 


21.07 


800.83 


28, 


.17 


28, 


.01 


28, 


.50 


28, 


.60 


29 


.20 


29, 


.41 


14 


83, 


,83023800 


-5, 


.32007800 


19.08 


22.96 


782.92 


28, 


.25 


28, 


.09 


28, 


.58 


28, 


.65 


29, 


.25 


29, 


.45 


15 


83, 


,83193200 


-5, 


.29617200 


27.37 


31.85 


764.63 


28, 


.42 


28, 


.25 


28 


.74 


28, 


.80 


29, 


.40 


29, 


.61 


18 


83, 


,83636500 


-5, 


.22322200 


2873.07 


4004.38 


698.35 


30 


.48 


30 


.31 


30 


.80 


30 


.94 


31 


.54 


31 


.75 


20 


83, 


,84798400 


-5, 


.46351500 


19.52 


23.76 


784.75 


28, 


.26 


28, 


.09 


28 


.58 


28, 


.66 


29, 


.26 


29, 


.47 


22 


83, 


,86342600 


-5, 


.44756500 


16.24 


20.93 


732.72 


28, 


.21 


28, 


.04 


28 


.53 


28, 


.63 


29, 


.23 


29, 


.44 


23 


83, 


,89975700 


-5, 


.35324100 


19.93 


23.70 


651.91 


28, 


.35 


28, 


.18 


28 


.67 


28, 


.74 


29, 


.34 


29, 


.55 


25 


83, 


,75458500 


-5, 


.39371600 


14.13 


18.14 


762.24 


28, 


.16 


27, 


.94 


28 


.37 


28, 


.56 


29, 


.15 


29 


.36 


32 


83, 


,80584286 


-5, 


.39145145 


14.85 


19.82 


346.70 


28, 


.52 


28, 


.30 


28 


.73 


28, 


.94 


29, 


.54 


29 


.74 


35 


83, 


,80742228 


-5, 


.39444168 


14.63 


22.59 


336.57 


28, 


.53 


28, 


.31 


28 


.74 


29, 


.01 


29, 


.60 


29, 


.81 


37 


83, 


,81955700 


-5, 


.31254300 


846.89 


1300.70 


780.42 


29, 


.92 


29, 


.71 


30 


.14 


30 


.40 


31 


.00 


31 


.21 



Table 11: X-ray parameters of the Class II candidates. This table is available in its entirety in electronic 
form; 



Seq. 


COUP 


log Nh 


kTl 


log EMI 


kT2 


log EM2 


log Lh.c 


log it,c 


Num. 


ID 


(cm~^) 


(keV) 


(cm~^) 


(keV) 


(cm~^) 


(ergs^i) 


(ergs"i) 


46 


12 


21.33± 0.18 


14.34± 11.68 


52.44± 0.05 






29.48 


29.64 


47 


13 


20.00± 2.12 


12.58± 15.00 


51.51± 0.17 






28.54 


28.71 


48 


21 


20.79± 2.06 


0.68± 0.11 


51.63± 0.17 


1.77± 0.47 


51.89± 0.12 


28.38 


29.04 


49 


29 


21.59± 0.03 


0.81± 0.06 


52.56± 0.13 


3.13± 0.32 


52.82± 0.03 


29.59 


30.04 


52 


37 


20.00± 2.34 


0.67± 0.33 


51.13± 0.63 






26.70 


28.10 


53 


40 


20.35± 2.00 


0.86± 0.11 


51.55± 0.15 


6.29± 4.63 


51.64± 0.08 


28.62 


28.99 


54 


44 


21.79± 0.32 


1.14± 0.35 


51.50± 0.29 






27.64 


28.44 



Table 12: X-ray parameters of the Class III candidates. This table is available in its entirety in electronic 
forrcL 



Seq. 


COUP 


log Nh 


kTl 


log EMI 


kT2 


log EM2 


log Lh,c 


log Lt,c 


Num. 


ID 


(cm~^) 


(keV) 


(cm~^) 


(keV) 


(cm~^) 


(ergs-i) 


(ergs-i) 


194 


6 


20.67± 0.43 


0.85± 0.12 


52.17± 0.11 


4.37± 0.95 


52.62± 0.04 


29.50 


29.84 


197 


7 


21.12± 0.12 


0.78± 0.09 


52.16± 0.12 


2.32± 0.27 


52.68± 0.04 


29.32 


29.79 


201 


9 


20.93± 0.04 


0.83± 0.01 


53.58± 0.02 


2.30± 0.08 


53.94± 0.01 


30.57 


31.08 


202 


10 


21.37± 0.03 


0.83± 0.02 


53.59± 0.03 


3.04± 0.11 


53.87± 0.01 


30.64 


31.08 


203 


14 


20.00± 1.21 


0.74± 0.13 


51.42± 0.21 






27.11 


28.39 


205 


20 


20.81± 0.14 


0.69± 6.24 


51.64± 1.00 


9.23± 15.00 


51.19± 0.21 


28.25 


28.82 
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Table 13: Upper Limits of the X-ray undetected Class II and III candidates 



Seq. 


ID 


RA (2000) 


Dec (2000) 


Cts 


Cts 


Exp. log L' 


log 


log LU 


log Ll 


c log LI 


c log Lie Class 


Num. 


H97 


(deg) 


(deg) 


hard 


tot. 


(ks) 


(erg s ^ 


(erg s ^ 


(erg s ^ 


(ergs- 


^ (ergs~ 


^ (ergs ^ 




51 


85 


83.69498179 


-5.35801543 


117 


130 


0.76 


28.83 


27.89 


29.10 


28.73 


28.54 


29.38 


II 


55 


108 


83.70820425 


-5.31230548 


83 


66 


0.45 


28.92 


27.99 


29.19 


28.66 


28.47 


29.31 


II 


66 


179 


83.74109068 


-5.47828377 


79 


237 


0.73 


28.67 


27.73 


28.94 


29.00 


28.81 


29.66 


II 


74 


219 


83.75862700 


-5.30628963 


27 


35 


0.77 


28.19 


27.26 


28.46 


28.15 


27.96 


28.81 


II 


98 


294 


83.77694276 


-5.51927167 


32 


39 


0.71 


28.29 


27.35 


28.56 


28.23 


28.04 


28.89 


II 


113 


402 


83.80437902 


-5.56763586 


57 


66 


0.66 


28.58 


27.64 


28.85 


28.50 


28.31 


29.15 


II 


123 


449 


83.81012490 


-5.55521463 


51 


500 


0.68 


28.51 


27.57 


28.78 


29.36 


29.17 


30.01 


II 


153 


682 


83.83838625 


-5.55481843 


47 


57 


0.71 


28.46 


27.52 


28.73 


28.40 


28.21 


29.05 


II 


166 


789 


83.85709585 


-5.49312187 


24 


29 


0.76 


28.14 


27.20 


28.41 


28.08 


27.89 


28.73 


II 


339 


729 


83.84592814 


-5.49483049 23.07 28.00 


0.72 


28.12 


27.21 


28.41 


27.92 


27.85 


28.77 


III 


210 


62 


83.67848933 


-5.42114559 38.07 55.00 


0.75 


28.31 


27.40 


28.60 


28.19 


28.12 


29.04 


III 



Table 14: Literature Class O/I counterparts. 
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